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Carrier density and temperature-dependent resistivity of graphene grown by chemical vapor depo- 
sition (CVD) is investigated. We observe in low mobility CVD graphene device a generic insulating 
behavior at low temperatures, and eventually a metallic behavior at high temperatures, manifest- 
ing a non-monotonic temperature dependent resistivity. This feature is strongly affected by carrier 
density modulation with the low-density samples exhibiting insulating-like temperature dependence 
upto higher temperatures than the corresponding high-density samples. To explain the tempera- 
ture and density dependence of the resistivity, we introduce thermal activation of charge carriers in 
electron-hole puddles induced by randomly distributed charged impurities. Our observed tempera- 
ture evolution of resistivity is then understood from the competition among thermal activation of 
charge carriers, temperature-dependent screening, and phonon scattering effects. Our experimental 
results are in good agreement with recent theories of graphene transport. 

PACS numbers: 72.80.Vp, 81.05.ue, 72.10.-d, 73.22.Pr 
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I. INTRODUCTION 

Graphene, a two dimensional honeycomb lattice of car- 
bon atoms composed of two identical interpenetrating tri- 
angular sub-lattices, has attracted research interest be- 
cause of its novel linear dispersion relation resulting in 
unique electronic properties such as room temperature 
intrinsic carrier mobility 1-4 . Exploiting such superior 
electronic properties combined with recently developed 
large scale growth of graphene by Chemical Vapor Depo- 
sition (CVD)£~— , has led to field effect devices integrated 
on a wafer scale showing promise for future electronic 
applications 8 . However, field effect performance of exfo- 
liated graphene is adversely affected by carrier scatter- 
ing by charged impurities in the surrounding media^^, 
surface phononsiir— on gate dielectric, and graphene 
acoustic phono n 11 ' 14 . Although graphene grown by CVD 
shows outstanding uniformity and mobility on a wafer 
scale, field effect transport measurement performed af- 
ter the transfer process, may be strongly affected by un- 
desirable impurity scattering^. Therefore, understand- 
ing transport on CVD-grown graphene influenced by im- 
purity scattering is an important issue for funamcntal 
physics and for further improvement of device perfor- 
mance. 

In this work we are presenting an experimental study 
of monolayer graphene transport properties, with the em- 
phasis on understanding the temperature dependent re- 
sistivity, complemented by a detailed comparison with 
theoretical calculations of graphene resistivity limited by 
charged impurity and phonon scattering. Such a com- 
bined experimental and theoretical study has not been 
undertaken for CVD graphene. Our work is closest 
in spirit to the early work by Tan et al. 16 in exfoli- 



ated graphene where, however, only the low-temperature 
density-dependent transport properties of graphene were 
studied in a combined experimental-theoretical investiga- 
tion to establish the dominant role of charged impurity 
scattering in determining the low-temperature graphene 
electrical conductivity. Early work on graphene trans- 
port properties 2,3 emphasized the weak temperature de- 
pendence of graphene resistivity upto the room temper- 
ature, presumably because of the weakness of the basic 
electron-phonon coupling in graphene. Such a weak tem- 
perature dependent conductivity is unheard of in elec- 
tronic materials where phonons invariably dominate the 
room-temperature resistivity. In graphene, the very weak 
phonon effects on the electrical conductivity have only 
recently been experimentally observed^ verifying earlier 
theoretical predictions^. Although there have been sev- 
eral earlier experimental works on the temperature de- 
pendence of transport properties in exfoliated graphene, 
our current work is the only work in the literature, to 
the best of our knowledge, which reports on the mea- 
sured temperature dependence of graphene resistivity in 
CVD graphene. Our work is also the only one in the 
literature combining experiment with detailed theory to 
shed light on the physical processes controlling the tem- 
perature dependence of graphene resistivity. One of our 
most important findings is that the low carrier den- 
sity (and low temperature) temperature-dependence in 
graphene is dominated to a large extent by carrier activa- 
tion across the potential fluctuations created by random 
charged impurities in the environment (which also lead 
to the inhomogeneous electron-hole puddle formation re- 
sponsible for the graphene minimum conductivity phe- 
nomenon). We compare our experimental results with a 
recent theory^ on the puddle induced temperature de- 
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pendence in graphene transport properties, obtaining 
excellent agreement. Such a comparison between theory 
and experiment for graphene transport at low tempera- 
ture (and densities) where puddles dominate transport 
properties is completely new in the literature. 

One other important indirect consequence of our work 
is the conclusion that graphene resistivity seems to man- 
ifest similar qualitative temperature and density depen- 
dence for similar mobility samples independent of how 
it is made, i.e. our CVD-grown samples exhibit a tem- 
perature dependent resistivity qualitatively similar to ex- 
foliation graphene samples of similar mobility. This es- 
tablishes beyond any doubt that the factor determining 
graphene transport properties (both density and tem- 
perature dependence) is the effective "zero-temperature" 
(i.e. very low-temperature) mobility which itself is de- 
termined by the charged impurity distribution in the en- 
vironment. Thus, improving graphene mobility for elec- 
tronic applications depends entirely on making graphene 
purer by removing the impurity content from the envi- 
ronment. Having CVD graphene helps only if it has low 
impurity content, which, according to our work, is no 
longer the case when CVD graphene is transferred. 

We investigate temperature (T) dependent resistivity, 
p{T), of graphene grown by CVD at various carrier den- 
sities (n g ). In low mobility samples, a strong insulat- 
ing behavior at low temperature and a metallic behav- 
ior at high temperature is observed, resulting in non- 
monotonicity. (We define insulating/metallic as dp/dT 
being negative/positive in the sample - this is purely an 
operational definition which is used extensively in the 
literature: all our samples are metallic in reality since a 
mobility can be defined.) This feature persists up to rel- 
atively high carrier density. On the other hand, high mo- 
bility samples show overall metallic behavior from fairly 
low to high carrier density. To understand this pecu- 
liar transport property of CVD graphene, we introduce 
the concept of thermal activation of charge carriers in 
electron-hole puddles generated by charged impurity 19 
which strongly affects the graphene mobility. Anoma- 
lous T and n g dependence of transport is well described 
theoretically by the competition among thermal activa- 
tion, screened impurity scattering, and phonon scatter- 
ing. Thus, puddles play a very important role in control- 
ling the temperature dependent CVD graphene trans- 
port. We compare our data with a recently developed 
detailed theory^. 

Temperature dependence of carrier transport proper- 
ties, e.g. resistivity, in electronic materials is often at- 
tributed to be arising from carrier-phonon interaction 
since phonons are thermally excited quantized- vibrations 
of the underlying lattice^. For example, the tempera- 
ture dependence of electrical resistivity in normal met- 
als like Al or Cu is essentially entirely controlled by 
the electron-phonon interaction strength in the metal. 
Such 'metallic' temperature dependence always involves 
an increasing resistivity with increasing temperature, i.e. 
dp/dT > 0, and often, the definition of a metal is taken 



to be synonymous with the dp/dT > behavior (or very 
weak temperature dependence of resistivity) with a finite 
(and presumably, not too large) zero-temperature extrap- 
olated resistivity. An insulator, on the other hand, has 
infinite (zero) resistivity (conductivity) at T = 0, and 
increasing T typically increases its conductivity since ac- 
tivated or variable-range hopping transport become al- 
lowed at finite temperature. An insulator or a localized 
system thus has dp/dT < with an extrapolated zero- 
temperature resistivity which is infinite (or very large). 
Thus, dp/dT < is often taken to be an operational 
definition of an insulator. This operational definition of 
an effective 'metal' (dp/dT > 0) versus an effective 'in- 
sulator' (dp/dT < 0) is general and powerful, and has 
long guided experimental and theoretical research in elec- 
tronic materials (metals, semiconductors, doped materi- 
als). However, we emphasize that 'metallic' and 'insu- 
lating' used in our paper are terminologies for the tem- 
perature dependence of resistivity. (dp/dT being nega- 
tive/positive does not necessarily indicate a T = true 
insulator /metal. ) 

Given the above very general consideration, it is there- 
fore always very interesting and intriguing when an oth- 
erwise well-defined metallic system exhibits insulating 
(i.e. dp/dT < 0) temperature dependence. It turns 
out that graphene is indeed such a rather unusual elec- 
tronic material which often manifests dp/dT < in 
some carrier density and temperature range while ex- 
hibiting the expected 'metallic' behavior dp/dT > in 
other regimes of density and temperature^ ^ 12 ' 14 . In 
fact, for some samples (and at some specific values of 
carrier density) even non-monotonic temperature depen- 
dence with resistivity decreasing with T first and then 
increasing with T is observed. In the current paper, we 
present a detailed experimental study of the temperature 
dependence of graphene transport properties comparing 
it with a quantitative theoretical analysis to show how 
and why such interesting temperature dependence could 
arise. The non-monotonicity in the graphene resistivity 
indicates competing physical processes, and indeed, there 
are several independent contributions to graphene trans- 
port which all depend non-trivially on density and tem- 
perature. These physical processes combine to give rise 
to the observed temperature dependence as we show in 
the current work. We emphasize that although some as- 
pects of this temperature dependence were earlier studied 
experimentally in the literature for exfoliated graphene 
on substrates and for suspended graphen o 11 ' 12 ' 20 ' 21 , ours 
is the first such temperature-dependent transport study 
for graphene grown by the CVD technique. Ours is also 
the first experimental study where the measurements are 
carefully compared with a detailed theoty. 

The competing physical mechanisms contributing to 
the temperature dependent graphene transport are 
phonons (which are known to be weak in graphene), 
screening, Fermi surface effects, and carrier activation 
across potential fluctuations associated with the density 
inhomogeneities or puddles in graphene which dominate 
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its low-density properties. Of these four mechanisms, 
the first two (phonons and screening) always produce 
metallic behavior whereas the last two (Fermi surface 
and activation) always produce 'insulating' behavior. In 
general, both insulating effects (i.e. Fermi surface ef- 
fects and activation across puddles) are stronger at low 
carrier densities whereas the two 'metallic' mechanisms 
are stronger at higher carrier densities. Thus, we expect 
graphene to preferentially exhibit metallic (insulating) 
behavior at higher (lower) carrier densities. Of course, at 
high enough temperature, phonons will eventually dom- 
inate, and all graphene samples should eventually ex- 
hibit metallic transport properties (i.e. dp/dT > 0), 
and therefore, at some densities we expect to see non- 
monotonicity in the graphene resistivity as a function 
of temperature. All these expectations are beautifully 
manifested in our experiments and we find reasonable 
quantitative agreement with theory. It is interesting and 
nontrivial, but eventually quantitatively understandable, 
that graphene shows a complex interplay of effective in- 
sulating and metallic behaviors superposed on each other 
as a function of density with the insulating (metallic) be- 
havior dominating low (high) densities without any real 
localization or gapped insulator physics being present in 
the system. Ours is the first complete study of this phe- 
nomenon in CVD graphene with direct comparison be- 
tween experiment and theory. 

The paper is organized as follows. In Sec. II we de- 
scribe the sample preparation and experimental results 
for resistivity. In Sec. Ill we provide the theoretical anal- 
ysis of measured non-monotonic resistivity. We conclude 
in Sec. IV with discussion. 



II. EXPERIMENTS 

Graphene devices in this study are synthesized by 
Inductively Coupled Plasma enhanced Chemical Va- 
por Deposition (ICP-CVD) on Cu substrate. Dur- 
ing growth process, the substrate is heated to 650°C 
within 10 min under ~ 1CP 7 torr, then treated with 
H2 plasma. After purging with Ar for a couple of 
minutes, C 2 H 2 is added (C 2 H 2 : Ar =1: 40) for 
graphene growth at the same temperature. For the 
graphene transfer, Graphene/Metal/Si0 2 /Si substrate 
was spin-coated with PMMA (Aldrich, 950 A4) and 
attached a pressure sensitive adhesive ultraviolet-tape. 
Peeling the tape against Si wafer physically separates 
Tape/PMMA/Graphene/Metal layer due to poor adhe- 
sion of metal film and Si0 2 . After etching of under- 
lying Ni/Cu by soaking in FeCl3 and cleaning in wa- 
ter, Tape/PMMA/Graphene layer was pressed onto the 
SiO 2 (300nm)/Si substrate with pre-patterned electrodes 
(Cr/Au of 10/50nm). The successive removal of tape 
and PMMA in methanol and acetone, respectively, leaves 
only graphene layer over pre-patterned marks. Graphene 
is etched to rectangular shape by 2 plasma as shown in 
the inset of Fig. 1(a). The distance between source and 



drain electrodes is 8/zm and the width of channel is 7/im. 
The T-dependent resistivity is recorded using 4-probe 
measurement technique for 4K < T < 300K under ~10- 
5 torr after 24 hour annealing at 150 °C. This annealing 
process before measurement gets rid of hole-doping, and 
moves the Dirac Point to near zero gate voltage (V g ). 

Single layer graphene was confirmed by Raman 
spectroscopy 2 ^ with single Lorentzian fit of 2D peak at 
2700 cm -1 in Fig. 1(a) and optical transmittance of 
97.7% at 550nm2£ in Fig. 1(b). The optical image and 
Raman mapping from our previous worki^ verified that 
transferred CVD-grown graphene has no sign of wrinkles 
or bi, tri-layers of graphene. According to Matthiessens 
rule, we assume the conductivity is the sum of contribu- 
tions from long-range scatterers and short range scatter- 
ers or residual resistivity. Hole mobility is then extracted 
using the following formula, 

fT(Vg)- 1 = - [p h C g (V g - h V g < Voirac, 

(1) 

where cr(V g ) is conductivity, /z/, is hole mobility, C g is 
capacitance of gate dielectric, Vdrag is the gate voltage 
at Dirac point, and u res ,h is residual conductivity. The 
measured hole mobility in a specific graphene sample in 
the same 6-inch wafer has a large variation by an order of 
magnitude. We speculate this large variation is coming 
from non-uniform transfer procedure which introduces 
local charge impurities differently. Two groups of de- 
vices with mobility of 1,000 ~ 3,000 cm 2 /Vs (A) and 
10,000 - 13,000 cm 2 /Vs (B) respectively were selected 
to understand the relationship between p(T) at various 
n g and the level of disorder (high in A, low in B^ilil. 
Fig. 1(c) and (d) show V g -dependent resistivity [p(V g )] 
of A and B devices, respectively. No marked difference is 
observed in Raman spectroscopy of these two disparate 
groups despite the order of magnitude variation in the 
mobility supporting introduction of non-uniform charge 
impurities during transfer. In addition, we observe a rel- 
atively weak D-peak, which indicates minimum level of 
structural defects in our samples. 

T-dependence of resistivity [p(T)] with reference to 
p(A.2K) at carrier densities near charge neutrality point 
is plotted in Fig. 2(a) and (b) for A and B respec- 
tively. AVg is equal to V g — Vbirac- The characteristic 
observed feature of carrier transport in our low mobil- 
ity sample A (Fig. 2(a)) is a strong insulating behavior 
(i.e., dp(T)/dT < 0) at low temperatures (T < 200K) 
for all densities (up to 7.2 x 10 n cra~ 2 ). For AV g > 4V 
(n g > 2.9 x 10 11 cm -2 ) a non-monotonicity in tempera- 
ture dependent resistivity develops (i.e. as T increases, 
p initially decreases and then increases forming a local 
minimum). The observed insulating behavior and non- 
monotonicity at high density in low mobility samples is 
not observed in high mobility samples B (Fig. 2(b)), 
where the measured resistivity shows a metallic temper- 
ature dependence at higher densities even though the ex- 
perimental conductivity manifests an insulating behavior 
at very low carrier density (up to n g ~ 2.0 x 10 n cm~ 2 ). 
Interestingly, negative dp(T)/dT at Dirac point is ~ 
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FIG. 1. (color online) (a) Raman spectroscopy on a device, 
inset: Optical microscope image of a device, (b) Light trans- 
mittance through CVD grown single layer graphene/Glass. 
(c) and (d) resistivity vs. gate voltage at various temperature 
for the device group of A and B, respectively. 



40% for both samples independent of sample mobility, 
which also has been reported in high mobility suspended 
graphen o 20 ' 24 . Recently, ballistic transport was invoked 
in order to explain the strong T-dependent p(T) at Dirac 
point as well as low n g regime 2 ^. However, this approach 
is not applicable in our low mobility devices which are ob- 
viously in the diffusive regime. The temperature depen- 
dence in our diffusive samples is presumably connected 
with extrinsic effects such as charged impurity scatter- 
ing, which is considered the dominant scattering process 
for transport at low n g . Note that only metallic behavior 
has been seen in graphene samples at high densities. 

Another important aspect of temperature dependent 
resistivity is that the initial insulating behavior at Dirac 
point is strongly affected by carrier density modulation. 
Starting from the charge neutrality point, the insulating 
behavior becomes non-monotonic near AV g ~ 4V~ (n g — 
2.9 x 10 n cm~ 2 ) for sample A due to the coexistence 
of low T insulating and high T metallic behavior (Fig. 
2(a)), but for sample B, a complete transition to metallic 



behavior occurs already around n g 
2(b)). 
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III. THEORETICAL ANALYSIS 

To understand T and n g dependence of different 
mobility devices, we consider the quantitative the- 
ory for carrier transport in the presence of electron- 
hole puddles due to the spatially fluctuating inhomoge- 
neous potential induced by random charged impurities. 
The formation of electron-hole puddles has been pro- 
posed theoretically 2 ^ and observed in scanning tunnel- 
ing spectroscopy 27 and scanning single electron transistor 
microscopy experiments 2 ®. The presence of electron-hole 
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FIG. 2. (color online) [p(T) - p{A.2K)]/ p{A.2K ) at varied 
AV 9 . (a) device A, (b) device B. 



puddles was recently introduced to explain the insulating 
behavior in metallic bilayer graphene^. It is crucial to 
include the role of electron-hole puddles to explain the 
measured insulating behavior in our monolayer graphene 
samples. We incorporate the puddles into the potential 
fluctuations which is assumed to be described by a statis- 
tical distribution function of a Gaussian form with a stan- 
dard deviation Thus, we calculate theoretically the 
carrier transport for our graphene samples by taking into 
account the highly inhomogeneous density and potential 
landscape. We also include the scattering by screened 
charged impurities, longitudinal acoustic phonon scatter- 
ing and surface optical phonon scattering. Then the total 
scattering time is calculated by 
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where r imp , r ap , and r sp are the transport scattering 
times of screened charged impurity 29 , longitudinal acous- 
tic phonon^, and surface optical phonon^, respectively. 
The scattering time due to charged impurity Ti mp is given 
by 
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where rij is the charged impurity density, £k = ftvfk 
is the carrier energy with the Fermi velocity vp and 
2D wave vector k, g(6*k,k') = [1 + cos(0)]/2 is a wave 
function form factor associated with the chiral matrix of 
graphene, 6*kk' is the scattering angle between k and k', 
Vi(q) = 2ne 2 / nq is the 2D Coulomb potential, and e(q) is 
the finite temperature dielectric screening function. The 
scattering time due to acoustic phonon mode r QC is given 
by 
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where f(e) = [eS £ ~^/ kBT + is the Fermi distribution 
function with the chemical potential /i, and Wkk' is the 
transition probability from the state with momentum k 
to k' state and given by 

WW = yj] |C(q)| 2 x [JV,*(e -e' + w q ) 
q 

+ (JV, + l)*(e-e'-w q )], (5) 

where C(q) is the matrix element for scattering by acous- 
tic phonon, w q = Up^q is the acoustic phonon frequency 
with v p h being the phonon velocity, and N q the phonon 
occupation number N q = l/(exp(/3w q ) — 1). The first 
(second) term is Eq. (0 corresponds to the absorp- 
tion (emission) of an acoustic phonon of wave vector 
q = k — k'. The matrix element C(q) for the deformation 
potential is given by 
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where D is the deformation potential coupling constant, 
p m is the graphene mass density, and A is the area of the 
sample. The scattering time due to the surface optical 
phonons is given by 



FIG. 3. (color online) p(T) at varied AV g . Scatters are 
experimental data and lines are calculation for (a) device A, 
(b) device B. 
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where the matrix elements of a surface optical phonon is 
given by 



[M Ql (q)} 2 = 



27re 2 



-2qd ^SOl 



1 



1 



1 



f 



(8) 



where d is the separation distance between graphene layer 
and substrate, eo (foo) is the static (high frequency) di- 
electric constant, and tosoi is the l-th kind of surface op- 
tical phonon frequency. Then by averaging over energy 
we have the conductivity as 
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where D(e) — 2e / Trlhvp) 2 is the density of states of 
graphene. 

We find that the screened charge impurity scattering 
dominates over phonon scattering at low T and strongly 
contributes to p(T) for all T. The acoustic phonon scat- 
tering gives rise to a linear T dependent resistivity and 
its contribution to p(T) is a quantitatively small even at 
room temperature. It is only important for high mobil- 
ity and high density samples because it is independent 
of n g and sample quality. It is known that the carriers 
in graphene are strongly coupled to the surface optical 
phonons of a polar substrate such as Si02- Consequently, 
the surface phonon scattering gives a large contribution 
to p(T) when the carrier temperature is high enough to 



absorb or emit phonons. We find that the surface phonon 
(^battering plays an important role for T > 200K and the 
scattering rate strongly depends on the distance between 
the 2D graphene layer and the substrate. 

We compare experimental and theoretical p(T) in Fig. 
3 incorporating all effects we discussed above. We use the 
charged impurity densities (m) and the standard devia- 
tion of potential fluctuation (s), (a) n, = 1.8 x 10 12 cm 2 
and s — 70 meV, (b) n, = 4.5 x 10 n cm 2 and s = 110 
meV. To get the charge impurity densities we first fit our 
data at low temperature and high density limits because 
in these limits all phonon scattering effects are exponc- 
tionally suppressed and scattering by charged impurity 
dominates. For acoustic phonon scattering we use the 
phonon velocity v P h = 2x 10 6 cm/s and the deformation 
potential coupling constant D = 19 eVJ^ For surface op- 
tical phonon scattering we use the two optical phonons in 
our calculation with uj sp \ — 63 meV and and lo SP 2 = 116 
meVfi 2 - The overall resistivity scale depends on the impu- 
rity density or mobility value, but the qualitative trends 
in p(T) arise from the basic aspects of the underlying 
scattering mechanisms. The rise in p with increasing 
T at high temperatures is a direct result of the ther- 
mal weakening of screening and phonon effects. The de- 
crease in p with temperature at low temperatures arises 
from thermal activation in the presence of electron-hole 
puddles 19 . The non-monotonicity of p(T) at the inter- 
mediate temperatures arises from the competition among 
activation, temperature dependent screening and phonon 
scattering. The activated insulating behavior competes 
with the metallic behavior induced by the temperature 
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dependent screening and phonon effects. When the po- 
tential fluctuation is smaller than the chemical potential 
the activation behavior is suppressed. As a result the to- 
tal conductivity manifests a metallic behavior. However, 
for large potential fluctuations the temperature depen- 
dence by activation overwhelms the temperature depen- 
dence due to all other mechanisms and the system shows 
insulating behavior. In our calculations we neglect An- 
derson localization. But we cannot rule out that in the 
lowest mobility sample localization may play a role in the 
insulating behavior. Detailed quantitative understanding 
in this regime would require a more sophisticated theory 
which includes higher order electron-electron interactions 
and disorder induced localization corrections, but our 
theoretical result o 13 ' 14 ' 19 ' 26 ' 29 give an excellent qualita- 
tive description of the data. We also measure the resis- 
tivity (not shown) for very low mobility sample (/z = 300 
cm 2 /Vs) showing a stronger insulating behavior at low 
temperatures than high mobility samples of Fig. 3(a) 
and (b). 

IV. DISCUSSION AND CONCLUSION 

In conclusion, temperature dependent transport of 
CVD grown graphene is investigated. For low mobility 
samples, an anomalous insulating behavior at low tem- 
peratures and non-monotonic at higher temperatures are 
observed, in contrast to the monotonic temperature de- 
pendent transport of high mobility exfoliation graphene. 
However, we note that a strong insulating behavior can 
be obsered at low mobility exfoliated graphene. Most 
important factor in observing an insulating temperature 
dependent transport is the strength of potential fluctu- 
ations which are proportional to the charged impurity 
density or inversely proportional to the mobility. If the 
sample mobility is low enough the insulating behavior is 
notable for even exfoliation graphene samples^. These 
are explained based on the transport theory incorpo- 
rating thermally activated transport of inhomogeneous 
electron-hole puddles. Our results imply that carrier 
transport in CVD graphene devices, which are vulnerable 
to the impurities or charged defects during the transfer 
process, is strongly affected by underlying environments, 



which gives rise to lower mobility in CVD graphene than 
the exfoliated graphene since both CVD and exfoliated 
graphene use the same substrate. To get higher mobility 
CVD graphene a careful transfer process is required in 
order to keep the background impurity density low. 

We have shown, through experimental measurements 
and detailed comparison with theoretical calculations, 
that the temperature dependence of graphene resistiv- 
ity can be understood as arising from a complex inter- 
play among the physical mechanisms of phonon scatter- 
ing, screening, Fermi surface averaging, and activation 
across potential barriers in inhomogeneous puddles. The 
first (last) two mechanisms lead to metallic (insulating) 
temperature dependence with resistivity increasing (de- 
creasing) with temperature, and therefore, the net effect 
could be either an insulating or a metallic temperature 
dependence or a non-monotonic combination of both, de- 
pending on the details of sample mobility, density, and 
the temperature regime. With increasing (decreasing) 
density and mobility, puddle effects are suppressed (en- 
hanced) , leading to stronger metallicity. With increasing 
temperature, phonon effects eventually dominate, so at 
the highest temperature (and not too low densities), the 
system always manifests metallic behavior with increas- 
ing resistivity with increasing temperature. At the Dirac 
point, activation in puddles dominate and the insulating 
behavior is generic. For intermediate temperatures and 
densities, the generic behavior of resistivity as a function 
of temperature is non-monotonic, as predicted in ref. [l9l . 
with the resistivity decreasing with increasing tempera- 
ture at low temperatures and then eventually increasing 
with temperature (except at the lowest densities where 
puddles are dominant even at room temperatures) at 
higher temperatures. The crossover temperature for the 
transition from insulating to metallic transport behavior 
increases with decreasing mobility and/or carrier density. 
All of these observed effects are completely consistent 
with the predictions of the theory developed recentl y 18 ' 19 . 
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